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Abstract 

Background: Similarly to the rest of the world, southern Africa's diverse chondrichthyan fauna is currently 
experiencing high fishing pressures from direct and non-direct fisheries to satisfy market demands for shark 
products such as fins and meat. In this study, the development of microsatellite markers through cross-species 
amplification of primer sets previously developed for closely related species is reported as an alternative approach 
to de novo marker development. This included the design of four microsatellite multiplex assays and their 
cross-species utility in genetic diversity analysis of southern African elasmobranchs. As this study forms part of a 
larger project on the development of genetic resources for commercially important and endemic southern African 
species, Mustelus mustelus was used as a candidate species for testing these multiplex assays in down-stream 
applications. 

Results: Thirty five microsatellite primer sets previously developed for five elasmobranch species were selected 
from literature for testing cross-species amplification in 16 elasmobranch species occurring in southern Africa. 
Cross-species amplification success rates ranged from 28.6%-71.4%. From the successfully amplified microsatellites, 
22 loci were selected and evaluated for levels of polymorphism, and four multiplex assays comprising of the 22 
microsatellites were successfully constructed, optimised and characterised in a panel of 87 Mustelus mustelus 
individuals. A total of 125 alleles were observed across all loci, with the number of alleles ranging from 3-12 alleles. 
Cross-species amplification of the four optimised multiplex assays was further tested on 1 1 commercially important 
and endemic southern African elasmobranch species. Percentage of polymorphism ranged from 31.8%-95.5% in 
these species with polymorphic information content decreasing exponentially with evolutionary distance from the 
source species. 

Conclusions: Cross-species amplification of the 35 microsatellites proved to be a time- and cost-effective approach 
to marker development in elasmobranchs and enabled the construction of four novel multiplex assays for characterising 
genetic diversity in a number of southern African elasmobranch species. This study successfully demonstrated the 
usefulness of these markers in down-stream applications such as genetic diversity assessment and species identification 
which could potentially aid in a more integrative, multidisciplinary approach to management and conservation of 
commercially important cosmopolitan and endemic elasmobranch species occurring in southern Africa. 

Keywords: Cross-species amplification, Microsatellites, Multiplex assays, Genetic diversity, Species identification, 
Conservation management 
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Background 

The subclass Elasmobranchii (sharks, skates and rays) 
comprises a diverse group of over 1000 species, and is 
representative of one of the most ancient extant verte- 
brate lineages [1]. Recently, pressures from direct and 
non-direct fisheries have resulted in the depletion of 
elasmobranch populations globally [2]. Decline in wild 
populations of elasmobranchs is further compounded by 
their life history traits that are more similar to those of 
mammals (e.g. low fecundity, late maturity and long ges- 
tation periods) than those of teleost fishes [3,4]. In com- 
parison, elasmobranchs may not respond well to the 
high fishing pressures. This trend has been particularly 
pronounced for sharks due to unregulated harvesting to 
support an increase in demand for shark products (e.g. 
fins, meat, liver oil, skin and cartilage). A drastic reduc- 
tion in population size (population bottleneck) can result 
in small populations experiencing the accumulating ef- 
fects of inbreeding leading to severe loss of genetic di- 
versity [5,6] . These trends have previously been reported 
for species such as the basking shark (Cetorhinus maximus) 
[7] and the narrownose smoothhound shark (Mustelus 
schmitti) [8]. Assessing genetic diversity and population 
structure of wild populations is therefore important for sus- 
tainable long-term management of the global shark fishery 
industry. 

Misidentification of shark species in fisheries opera- 
tions is also a widespread concern [9-12], and molecular 
individual identification methods have been developed 
to alleviate this problem [4,13-17]. To integrate genetic 
knowledge with fisheries management, it is imperative 
for shark fisheries to report shark landings by species in- 
stead of lumping them into species- or family groups 
(e.g. houndsharks, carcharinids, hammerheads etc.). This 
stems from the difficulties involved with unambiguously 
identifying species within and across families [13,18] e.g., 
carcharinids (Carcharhinus brachyurus, C. obscurus and 
C. plumbeus) and houndsharks (Mustelus mustelus, M. 
palumbes and Galeorhinus galeus) due to a high degree 
of conserved interspecific morphology [14]. Neglecting 
to report shark landings by species overlooks important 
differences in species susceptibility and population vul- 
nerability to exploitation [15], and that in turn has im- 
portant implications for species-specific conservation, 
management and trade monitoring programmes [19]. 

The general lack of molecular genetic markers (e.g. 
microsatellites) for many elasmobranch species impedes 
population and conservation genetic studies in that these 
markers can provide valuable information relating to 
population dynamics (spatial and temporal genetic vari- 
ation) of individual species. Microsatellites are highly poly- 
morphic due to their high mutation rate (between 10~ 3 
and 10" 4 mutations per gamete per generation) resulting 
in extensive length polymorphism [20,21]. This makes 



microsatellite markers one of the most powerful molecu- 
lar genetic tools with a remarkable array of applications 
ranging from genetic diversity [22,23] and population 
structure inference [24,25] to discerning genetic mating 
systems [26,27] and the identification of species [28-30] . 

Because the de novo development of microsatellites is 
challenging due to notoriously low rates of polymorphism 
in elasmobranchs [31], the development of microsatellite 
markers through cross-species amplification is the most 
effective alternative approach to de novo development of 
microsatellites and has recently also been reported in 
sharks [32]. Microsatellite cross-species amplification re- 
lies on the presence of conserved microsatellite flanking 
sequences [33], which in some organisms markedly dem- 
onstrate a high degree of conservation following millions 
of years of divergent evolution (e.g., 250 million years in 
sharks [28] and 470 million years in fish [34]). The success 
rate of microsatellite cross-species amplification has dir- 
ectly been correlated to the evolutionary distance between 
the source species and the target species [33,35]. 

The underrepresentation of endemic taxa in many 
cross-species amplification studies is unfortunate as en- 
demics should be of great interest for conservation of 
biodiversity on a regional scale. Southern Africa has one 
of the most diverse chondrichthyan faunas in the world, 
consisting of some 181 species in 44 families of which 
34 species are endemic to southern Africa [36,37]. Grow- 
ing concerns regarding the sustainability of the southern 
African shark fishery, stemming from the local declines of 
cosmopolitan and endemic species, lead to stricter regula- 
tions being imposed so as to avert the collapse of natural 
populations [38-40]. Accordingly, we report here the de- 
velopment of microsatellite markers through cross-species 
amplification of species-specific primers from closely re- 
lated species. This included the design and optimisation of 
four microsatellite multiplex assays and their cross-species 
utility in genetic diversity analysis of 11 southern African 
elasmobranch species. 

Results and discussion 

Cross-species amplification 

Development of microsatellite loci through cross-species 
amplification proved useful in establishing genetic markers 
for shark species that are commercially important and 
those (typically endemics) that are indirectly affected by 
fisheries' operations. Amplification of the 35 microsatellites 
in 50 individuals from 16 different elasmobranch species 
(1-4 individuals per species) proved to be effective (Table 1). 
Cross-species amplification success rates or the percent- 
age of microsatellites that amplified successfully ranged 
from 60.00%-71.40% in the Triakidae and Carcharhinidae 
families and from 28.57%-48.57% in the Scyliorhinidae, 
Sphyrnidae and Rajidae families (Figure 1). The higher 
success rates in the Triakidae and Carcharhinidae is 
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Table 1 Cross-species amplification of the 35 microsatellites among 16 elasmobranch species of southern Africa 
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-, no visible band or faint bands with insufficient band intensity for scoring alleles were observed; +, solid bands with sufficient intensity for scoring alleles were 
detected; ++, solid bands with artefacts were produced but with at least one band of expected allele size. Mustelus mustelus (MM), Mustelus palumbes (MP), 
Galeorhinus galeus (GG), Scylliogaleus quecketti (SQ), Carcharhinus brachyurus (CB), Carcharhinus limbatus (CL), Carcharhinus obscurus (CO), Carcharhinus plumbeus 
(CP), Haploblepharus pictus (HP), Haploblepharus edwardsii (HE), Poroderma africanum (PA), Poroderma pantherinum (PP) Sphyrna lewini (SL), Sphyrna zygaena (SZ), 
Raja straeleni (RS) and Raja alba (RA). 



expected as most microsatellites tested in this study were 
originally developed for species within the Triakidae fam- 
ily. Overall, the microsatellites showed less successful 
cross-species amplification to the taxa more divergent 
from the source species. Notably none of the individuals 



showed PCR amplification at any of the six Scyliorhinus 
canacula microsatellites. This may in part be attributed to 
S. canacula being more distantly related to the study spe- 
cies. The mean genetic distance between the taxa was 
21.4 ± 1.7% (mean ± SD) (G. galeus as source species; 
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Figure 1 Success rates of 35 microsatellite loci across five families of southern African elasmobranch species. 



Figure 2) and 18.7 ± 1.5% (M canis as source species; 
Figure 3). Haploblepharus pictus could not be repre- 
sented in the distance plot due to the lack of genetic in- 
formation available in GenBank and Global Cestode 
Database: Elasmobranchs Specimens. 

Results of cross-species amplification performance of 
G. galeus microsatellites exhibited a logarithmic regres- 
sion function (Figure 2), albeit non- significant (slope 
within the 95% CI for no difference from zero), that may 
explain the general trend of negative correlation between 
cross-species amplification performance and genetic di- 
vergence seen across taxa [30,33,35,41]. Cross-species 
amplification of microsatellite markers from source to 
target species is generally negatively correlated with evo- 
lutionary divergence [30,42-45]. A similar trend was not 
observed for the M. canis microsatellites (Figure 3). This 
can be due to different life history traits {i.e. mating sys- 
tem and generation time) and genome size (C value) be- 
tween the source and target species, which have been 
previously found to have significant negative effects on 
cross-species amplification success [33]. However, apart 
from the source-target species evolutionary distance 
other factors, such as mutations in microsatellite flank- 
ing sequences, may affect the success rate of cross- 
species amplification. Since microsatellites are usually 
found in non-coding regions where the substitution rate 
is higher than in coding regions [46], these microsatellite 
flanking sequences which serve as regions for PCR pri- 
mer design and binding sites are prone to mutations 
[35]. Mutations (indels) in these regions may therefore 
result in null alleles and in turn affect the patterns of 
cross-species amplification as demonstrated in birds [47] 
and salmonids [48]. 



Additionally, M. canis microsatellite loci were isolated 
from an enriched genomic library [49] whereas for G. 
galeus the microsatellites were developed by a high- 
throughput sequencing approach (Roche 454 pyrosequenc- 
ing) [50]. Based on the observed data it is hypothesised 
that these different approaches may have influenced the 
cross-species performance possibly due to the different 
resolving power of each approach for capturing micro- 
satellites distributed across different parts of the gen- 
ome. Castoe et al. [51] argues that enrichment-based 
approaches commonly use a few specific repeated mo- 
tifs, which are largely selected without prior knowledge 
of their abundance in the genome and therefore could 
introduce potential bias in genome representativeness. 
In contrast, microsatellite identification from randomly 
sequenced genomic regions {e.g. Roche 454 NGS) allows 
for an unbiased assessment of all types of microsatellite 
loci present in a genome [51]. 

Cross-species amplification of orthologous microsatel- 
lites, due to the presence of conserved microsatellite 
flanking sequences [33], can persist over millions of 
years following divergent evolution as previously shown 
in sharks (250 million years [28]) and in fish (470 million 
years [34]). This indicates lower mutation rates within 
microsatellite flanking regions in aquatic organisms [30,34]. 
The reported cross-amplified microsatellite markers will 
thus contribute to establishing a molecular genetic marker 
repository for each of the Southern Africa elasmobranchs 
species included in this study. Future research efforts may 
be dedicated to generating microsatellite primers that have 
a high cross-species utility {e.g. [52]) as well as in silico min- 
ing of polymorphic microsatellite markers from expressed 
sequence tag data [53]. 
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Figure 2 Cross-species amplification performance of Galeorhinus galeus microsatellites in 15 of the 16 elasmobranch species, and 
genetic divergence between G. galeus and the target species based on ND2 sequences. 



Multiplex assay characterisation 

Twenty- two microsatellite loci that successfully cross - 
amplified across the study species and showed poly- 
morphism in initial screening tests were used to develop 
four multiplex assays (MPS) comprised of at least five 
microsatellites each. These multiplex assays were charac- 
terised in a panel of 87 common smoothhound sharks 
(Mustelus mustelus) (Table 2). All 22 microsatellite loci 
were polymorphic across the multiplexes. With the 



exception of one locus, McaB22, all the microsatellite 
loci were found to deviate significantly from Hardy- 
Weinberg equilibrium most likely due to Wahlund effect 
as samples were pooled from diverse geographical loca- 
tions for analysis. MicroChecker detected no significant 
genotyping errors but indicated that null alleles were 
present at two loci (Mh9 and Gg7). Slatkins exact test 
for neutrality indicated that two loci (McaB22 and Gg3) 
were candidates for being subjected to selection. 
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Figure 3 Cross-species amplification performance of Mustelus canis microsatellites in 15 of the 16 elasmobranch species, and genetic 
divergence between M. canis and the target species based on ND2 sequences. 
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Table 2 Characterisation of four multiplex assays for Mustelus mustelus based on 87 individuals from southern Africa 
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0.455 


-0.733** 


-0.263 


0.488 


Gg23 


(AQ n 


0.2 


VIC 


258-278 


6 


2.84 


1.000 


0.651 


0.582 


-0.540** 


-0.214 


0.562 


MPS4 (mean) 










4.4 


2.08 


0.788 


0.461 


0.3816 


-0.589 


-0.211 


0.615 


Overall (mean) 










5.7 


2.4 


0.691 


0.512 


0.444 


0.010 


-0.139 


0.561 



Primer concentration in the final reaction as uM/primer ([P]); Number of alleles per locus {A N ); effective number of alleles {A E ); observed heterozygosity (H G ); 
expected heterozygosity (H £ ); polymorphic information content (P/C); inbreeding coefficient (F /s ) with statistically significant deviations from Hardy-Weinberg 
expectations indicated by *{P< 0.01) and **{P< 0.001); null allele frequency {Fr NULL ) with indicating the presence of null alleles at statistical significance at the 5% 
nominal level and Ewans-Watterson probability (P E . W ). Mean values for each multiplex assay and overall are indicated in bold. 



Multiplex assay cross-species amplification and efficiency 
in species identification 

Cross -species amplification of the four multiplex assays 
was tested for 11 other southern African shark species 
(Table 3). The number of alleles observed was highest in 
G galeus and M. palumbes, varying from 1 to 7, while 
the percentage of polymorphism (PP) for each marker 
ranged from 31.896-95.5%. The polymorphic information 
content {PIC) decreased exponentially with evolutionary 
distance from the source species (Table 3) and the four 
multiplex assays showed the highest PIC in M. mustelus, 
M. palumbes and G. galeus. 

The mean genetic diversity estimates for each species 
in terms of number of alleles (A N ), effective number 
of alleles (A E ), observed heterozygosity (H Q ), expected 
heterozygosity (H E ) and PIC are shown in Figure 4. In 



group 1, the mean H E , A E and PIC varied from relatively 
low in C. brachyurus (mean H E = 0.230; mean A E = 1.4; 
mean PIC = 0.181) to relatively high in M. palumbes 
(mean H E = 0.653; mean A E = 3.3; mean PIC = 0.606). 
Group 2 exhibited similar patterns of genetic diversity 
that varied from moderate in S. zygaena (mean H E = 
0.593; mean A E - 3.2; mean PIC = 0.554) to relatively 
high in P. pantherinum (mean H E = 0.662; mean A E = 
3.4; mean PIC = 0.603). For group 3 with n = 4, the mean 
H E) A E and PIC ranged from relatively low in C. plum- 
beus (mean H E = 0.249; mean A E = 1.5; mean PIC - 
0.193) to relatively high in C. obscurus (mean H E = 0.429; 
mean A E = 2.1; mean PIC = 0.367). 

The genotypic distribution of the study species is 
depicted in Figure 5. Most of the study species could be 
differentiated on PCI and PC2 of the PCoA plot as can 
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Table 3 Multiplex transferability results of a total of 22 microsatellite loci showing the number of alleles per locus for 
the 1 1 elasmobranch species tested 



Locus 


MP 


GG 


CB 


CL 


CO 


CP 


HP 


HE 


SL 


SZ 


PP 




(n = 8) 


(n = 8) 


(n = 8) 


(n = 4) 


(n = 4) 


fn = 4J 


(n =8) 


(n = 4) 


fn = 5J 


fn = 5J 


fn = 5J 


MPS1 
























Mhl 


4 


3 


1 


2 


1 


1 


1 


1 




2 


3 


Mh2 


- 


5 


* 


1 


1 


-X 


- 


- 


? 


1 


3 


Mh9 


4 


4 




- 


* 




- 


- 


4 


5 


3 


Mh25 


5 


5 


2 


3 


2 


3 


1 


3 


4 


6 


3 


Mca25 


3 


1 


1 


3 


2 


1 


2 


1 


3 


4 


3 


McaB39 


3 


3 


2 


2 


1 


2 


2 


* 


2 


3 


3 


MPS2 
























McaB5 


3 


1 


2 


4 


2 


3 


1 


1 


5 


5 


7 


McaB6 


4 


4 




4 


2 


1 


- 


1 


2 


3 


5 


McoB22 


2 


1 


1 


4 


2 


4 


2 


4 


4 


7 


7 


McoB27 


2 


2 


2 


- 


1 


-X- 


- 


1 


- 


2 


- 


Mca33 


4 


2 


2 


4 


2 


2 


2 


3 


4 


7 


6 


McoB37 


3 


5 


1 


3 


1 


1 


- 


- 


4 


6 


7 


MPS3 
























Gg2 


5 


4 


1 


- 


1 


1 


2 


- 


- 


2 


- 


Gg3 


3 


2 


1 


1 


2 


-X- 


2 


1 


- 


2 


2 


Gg7 


4 


1 


1 


- 


1 


1 


2 


- 


- 


1 


1 


Ggll 


6 


4 


1 


- 


2 


2 






- 


1 


3 


Ggl2 


4 


5 


- 


- 


1 


■X 


1 


1 


2 


2 


1 


MPS4 
























Gg15 


7 


5 


1 


4 


1 


4 


4 


2 


5 


3 


5 


Ggl7 


4 


4 


1 


4 


1 


2 


2 


2 


2 


3 


2 


Gg18 


6 


3 


2 


3 


1 


3 


3 


3 


7 


6 


4 


Gg22 


6 


5 


2 


3 


2 


2 


2 


7 


3 


2 


3 


Gg23 


4 


2 


2 


1 


2 


1 


3 


2 


6 


4 


3 


Total P loci 


21 


18 


8 


13 


10 


10 


12 


7 


16 


19 


18 


PP 


95.5 


81.8 


36.4 


59.1 


45.5 


45.5 


54.5 


31.8 


72.7 


86.3 


81.8 



n - number of individuals tested; *Failed to amplify but showed successful transferability initially (see Table 1); - No amplification; allele numbers in bold indicate 
loci that failed to cross-amplify according to Table 1; P - polymorphic and PP - percentage of polymorphism. For species abbreviations refer to Table 1. 



be seen from individuals of each respective species clus- 
tering together. Individuals of the catshark species 
(H. edwardsii, P. pantherinum and H. pictus) however, 
were dispersed across quadrant 3 and 4. The PCoA also 
revealed that one of the M. mustelus individuals was 
misidentified as G. galeus. The identity of this particular 
individual was subsequently confirmed using the genetic 
identification method developed specifically for smooth- 
hound sharks [16]. Briefly, this method involves using 
four primers (1 universal forward primer and 3 species- 
specific reverse primers) for the mitochondrial gene, 
nicotinamide adenine dehydrogenase subunit 2 (ND2), 
in a multiplex PCR reaction. The reverse primers amp- 
lify a fragment of different length for each species (M. 



aster ias, 564 bp; M. mustelus, 392 bp; G. galeus, 671 bp) 
and can therefore be utilised for distinguishing species 
based on fragment size. 

In this study, the potential use of microsatellite loci 
in species identification was successfully demonstrated 
using shared microsatellite loci between species. The 
polymorphic information of these microsatellite loci was 
characterised by low genetic variation as previously pro- 
posed for elasmobranchs [31]. The genotypic distribu- 
tion of the study species could also be differentiated 
based on PCoA analysis. Markedly, the lack of differenti- 
ation between the catshark species (H. edwardsii and 
H. pictus) on the PCoA plot may be explained by the 
misidentification of the Haploblepharus species that is a 
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common occurrence in the catshark family [54]. To 
further investigate whether the lack of differentiation de- 
tected with microsatellites was indeed due to misidentifi- 
cation, the cytochrome b (Cyt b) and ND2 genes were 
applied for species comparisons [54]. For both genes, se- 
quence analyses revealed individuals with cryptic identi- 
fication, suggesting that in the case of the catsharks, the 
microsatellites optimised in the current study were not 
successful in discriminating between the Haploblepharus 
species but could in the future aid in the identification 
of cryptic speciation within the catshark family. 

In South Africa, the aforementioned misidentification 
issue is prominent in fishing operations (particularly in 
longline and trawl fisheries) where there is a high rate 
of incidental capture of non-target shark species [9-12]. 
This hinders the collection of reliable data on shark 
catch and trade on a species-specific basis making ro- 
bust stock assessments and identification of overfished 
and potentially threatened species nearly impossible in 
most situations [4]. This was apparent in a study by 
Attwood et al [12], which assessed bycatch in South 
Africa's inshore trawl fishery based on observer records. 



In the aforementioned study, certain taxa were difficult to 
identify, and therefore taxonomic groups were lumped in 
species groups {e.g. Raja spp., Mustelus spp. and Scyliorhi- 
nidae), even though every attempt was made to analyse 
data at species level. Therefore, the molecular genetic 
makers developed by the current study may facilitate in 
obtaining species-specific catch data for stock assessment, 
characterising genetic diversity and delineating population 
genetic structure. This in turn will contribute to the im- 
plementation of future conservation and management 
plans on a species-specific level in southern Africa. 

Conclusions 

Cross-species amplification of available microsatellite 
loci to target species has proven to be more time- and 
cost-effective in comparison to the de novo development 
approach and permitted the cross-amplification of 22 
markers across 12 elasmobranch species. Cross-species 
amplification of the four multiplex assays developed in 
the current study highlighted the usefulness of microsa- 
tellites for characterising genetic diversity and potentially 
also species identification of a number of commercially 



• M. mustelus ■ M. palimibes a G. galeits * C. braclnwus 

♦ C. limbatus ■ C. obscwus * C. pi urn bens • H. pi cms 
H. edn ardsii ■ P. pantheriinmi S. lewini • S. zygaena 

Figure 5 Principle coordinates analysis (PCoA) of study species based on 12 shared amplified microsatellite loci between species. 

Arrows depict misidentified/mislabelled individuals. 
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important and endemic elasmobranch species. The mo- 
lecular genetic markers developed in this study and their 
usefulness in down-stream applications could therefore 
aid in a more integrative, multidisciplinary approach to 
conservation management of elasmobranchs in southern 
Africa. 

Methods 

Ethics statement 

The collection of specimens from various shark species 
used in this study complied with the Convention on Bio- 
logical Diversity (http://www.cbd.int/convention/) and 
the Convention on the Trade in Endangered Species of 
Wild Fauna and Flora (http://www.cites.org/). All per- 
mits to collect finclip or muscle tissue for research pur- 
poses were granted by the Department of Agriculture, 
Forestry and Fisheries (Republic of South Africa). 

Study species and DNA extraction 

Sixteen elasmobranch species occurring in southern 
African waters belonging to five families within two orders 
were selected for cross-species amplification (Additional 
file 1). Where possible, specimens were collected from at 
least two sampling locations to better capture allelic diver- 
sity present within populations of each respective species. 
However, due to opportunistic sampling for a majority 
of the study species, samples were obtained from only 
a single location (Additional file 1). Muscle tissue or 
finclips were preserved in 99% ethanol and stored at 
room temperature until further use. Total genomic DNA 
was isolated using the standard cetyltrimethylammonium 
bromide (CTAB) method of Saghai-Maroof et al. [55]. 
The extracted DNA was quantified using a NanoDrop 
ND-1000 spectrophotometer v.3.0.1 {NanoDrop"), For test- 
ing cross-species amplification, each DNA sample was ad- 
justed to a working concentration of 50 ng/ul and stored at 
-20°C prior to polymerase chain reaction (PCR) analysis. 

Microsatellite primer transfer 

A total of 35 microsatellite markers previously developed 
in five elasmobranch species (Raja pulchra [44], Scylior- 
hinus canacula [45], M. canis [49], G. galeus [50] 
and M. henlei [56,57]) were selected for testing cross- 
species amplification. Primer sequences and annealing 
temperature (T A ) of each primer set optimised for 
each respective source species are shown in Table 
S2 (Additional file 2). Polymerase chain reaction condi- 
tions optimised for the majority of the source species 
[G. galeus, M. henlei, R. pulchra and S. canacula] were 
applied for cross-species microsatellite examinations 
except for the M. canis ([56] PCR protocol) and some 
R. pulchra (modified [57] PCR protocol) primer sets. 
Polymerase chain reactions for all individuals were exe- 
cuted in a GeneAmp® PCR System 2700. 



The PCR amplicons were visualised on a 2% agarose 
gel stained with ethidium bromide together with nega- 
tive controls and Promega 100 bp molecular size ladder 
for preliminary size determination. Success or failure of 
PCR amplification in cross-species trials was determined 
simply on the basis of whether band intensity was suffi- 
cient to score alleles. In most instances, less stringent 
PCR conditions were not employed in the cross-species 
assays so as to minimise the risk of amplification of non- 
orthologous loci in the target species. The number of 
markers that showed amplification success at all or a 
percentage of individuals in the target species ("+/++" in 
Table 1) were counted as an index to measure the cross - 
species microsatellite amplification performance. 

Multiplex design and optimisation 

As this study forms part of a larger project on the devel- 
opment of genetic resources for commercially important 
and endemic species of southern Africa, Mustelus mus- 
telus was used as a candidate species for testing of the 
four novel multiplex assays. Levels of polymorphism were 
initially assessed at all the successfully cross-amplified 
microsatellite loci in a panel of eight individuals of 
M. mustelus. Amplicons were subjected to electrophoresis 
for two hours at 150 volts on a 12% polyacrylamide gel to 
detect size variants. Microsatellites were considered to be 
polymorphic when two bands were distinguishable in a 
single individual {i.e. heterozygote) and/or clear size differ- 
ences were detected between different individuals. 

Twenty-two polymorphic microsatellite loci were se- 
lected, and primers fluorescently labelled and optimised 
in four multiplex assays (5-6 loci per MPS) using a 
strategy outlined by Guichoux et al. [58] with one of the 
following dyes: FAM, VIC, PET, or NED. The use of dif- 
ferent dyes was to facilitate co-amplification of multiple 
microsatellite markers in a single reaction for cost- and 
time-efficient genotyping (Multiplex PCR). 

After optimisation of the newly designed MPS (MPS1, 
MPS2, MPS3 and MPS4), a panel of 87 M. mustelus indi- 
viduals from across the distribution range in southern 
Africa was genotyped for marker characterisation pur- 
poses. The multiplex assays were then tested on 11 add- 
itional species to show their overall application in genetic 
diversity and population structure analysis. Finally, a total 
of 12 microsatellite loci that were successfully genotyped 
across the study species (Mhl, Mh25, Mca25, McaB39, 
McaBS, McaB22, Mca33, GglS, Ggl7, Ggl8, Gg22 and 
Gg23) were selected to demonstrate the potential use of 
microsatellite loci in species identification. 

The percentage of polymorphism (PP) was calculated 
using the formula: 

N P 

PP = — X 100 

N T 
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where N P is the total number of polymorphic loci and 
N T is the total number of loci multiplied by 100. 

For the multiplex reaction, the Qiagen Multiplex PCR kit 
was used and PCR conducted according to the manufac- 
turers instructions except for varying T A , 59°C for MPS1, 
MPS3 and MPS4; and 56°C for MPS2. For subsequent ana- 
lysis on an ABI 3730XL DNA Analyzer, PCR products 
were diluted in distilled water and fragment analysis per- 
formed together with the LIZ600 internal size standard. 
Individual genotypes were scored based on fragment size 
via Peak Scanner® software v.l (Life Technologies). Auto 
Bin v.0.9 macro for Excel (http://www.bordeaux-aquitaine. 
inra.fr/biogeco/Ressources/Logiciels/Autobin; see [58]) was 
used to detect discreet size variants where allele binning of 
genotype data obtained from Peak Scanner® software v.l 
was based on raw size. 

Genetic diversity analysis 

MicroChecker v.2.2.3 [59] was used to evaluate the pres- 
ence of genotypic errors caused by allele dropout, stut- 
tering and null alleles. Null allele frequencies {Fr NULL ) 
were calculated using the Brookfield 1 estimator imple- 
mented in this program. Locus-specific fixation index 
(F IS ) and over all loci was estimated to measure departure 
from Hardy- Weinberg equilibrium using the exact prob- 
ability test (20 batches, Dememorization; 10000 and 5000 
iterations) using Genepop v.4.0 [60]. Linkage disequilib- 
rium between all pairs of loci was calculated using an 
exact test implemented also in Genepop. Slatkins exact 
test (1000 permutations) for neutrality, based on Ewens- 
Watterson sampling theory [61] was used to detect loci 
under selection as implemented in Arlequin v.3.5.1.2 [62]. 
The number (A N ) of alleles at each microsatellite locus, as 
well as the effective number of alleles A E : 



i/& 2 



i=l 



where p t is the frequency of the £ allele and n is the 
number of alleles was calculated using the GenAlEx 
v.6.5 program [63]. The proportion of individual samples 
that were heterozygous [direct count heterozygosity 
(H Q ) and expected under Hardy- Weinberg equilibrium 
(H E )] was calculated using MsatTools [64]. MsatTools 
was also used to calculate the polymorphic information 
content (PIC) of each marker according to the following 
equation in [65]: 

n-l n 



PIC = 1-^ Pi 2 -2 ^_^Pi 2 Pj 2 where p t and pj are 

_i=l j=i+l 

the frequency of the i th and / h allele respectively and 
n is the number of alleles. 

Direct comparison of genetic diversity estimates (H E} A E 
and PIC) across the 11 species was not plausible due to 



the different sample sizes that were used. Species were 
therefore grouped into three groups according to sample 
size: (1) M. mustelus, M. palumbes, G. galeus, C bra- 
chyurus and H. pictus (n = 8); (2) P. pantherinum, S. lewini 
and S. zygaena (n = 5) and (3) C. obscurus, C limbatus, C 
plumbeus and H. edwardsii (n = 4>). The potential use of 
microsatellite data for species-assignment was assessed 
through principle coordinate analysis (PCoA) in GenAlEx 
v.6.5 using genetic distances between individuals. 

To evaluate cross-species amplification performance, 
DNA sequences derived from the mitochondrial ND2 
gene (1044 bp) of each species were downloaded from 
GenBank and Global Cestode Database: Elasmobranchs 
Specimens (http://elasmobranchs.tapewormdb.uconn.edu) 
(Additional file 3). The genetic distance of the study taxa 
was estimated using the Kimura 2-parameter model with 
the rate variation among sites modelled with a gamma dis- 
tribution (shape parameter = 5) implemented in MEGA 
v.5 [66]. 

Additional files 



Additional file 1: The 16 elasmobranch species of southern Africa 
selected for cross-species amplification, including family, species, 
distribution and sampling locations. 

Additional file 2: The 35 putative microsatellite markers developed 
from five closely related species for cross-species amplification in the 
study taxa, including the primers sequence, microsatellite repeat motif, 
annealing temperature (T A ) and GenBank accession numbers. 

Additional file 3: The ND2 sequence information of the study taxa 
used to estimate the genetic distance to evaluate cross-species 
performance, including ID Verified, availability of images (yes or no) 
which are available in the on-line host specimen database (http:// 
elasmobranchs.tapewormdb.uconn.edu) and GenBank accession 
numbers. 
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